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Abstract 

We study in a quark model the contradiction between factorization and 
duality found in nonleptonic decays at next to leading order in 1/N C , concen- 
trating on quark exchange mechanism. The contradiction originates in the fact 
that the standard factorization assumption approximates the asymptotic final 
states by a non-orthogonal set of states, thus leading to an overcounting of the 
decay probability. We consider a system with two heavy quarks treated as clas- 
sical color sources with constant velocity, and two mass-degenerate antiquarks. 
Exploiting permutation symmetry in an adiabatic approximation, we find that 
final state interaction restores duality. Three 0(1/N C ) effects are exhibited: i) 
a proper treatment of orthogonality yields a global correction 1/N C — > 1/2N C 
within a generalized factorization a la BSW, (such a factor was present in an 
Ansatz by Shifman) , ii) the distorsion of the meson wave functions at the time 
of the weak decay, iii) relative phases generated by the later evolution. The 
latter effect becomes dominant for light antiquarks or for a small velocity of 
the final mesons, and may thoroughly modify the factorization picture. For 
exclusive decay it may interchange the role of class I and class II final channels, 
and for semi-inclusive decay it may lead to an equal sharing of the probability 
between the two sets of final states. In the heavy antiquark and large velocity 
limit, the replacement 1/N C — * 1/2N C is the dominant correction. 
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1 Introduction. 



While some progress has been performed during these last years in our understanding 
of semileptonic and leptonic decay mechanisms, our understanding of nonleptonic 
decays is still semi-quantitative. Not to speak of the AI = 1/2 mystery, the non- 
leptonic decays of D and B mesons are most often studied with the help of the 
standard factorization assumption the theoretical basis for which exists only in 
the N c — > oo limit, or of the generalized factorization assumption a la BSW H which 
is a phenomenological Ansatz. 

The nonleptonic decay channels of heavy meson are an important issue and 
will grow even more so, since they provide the cannels in which the CP asymmetries 
will be looked for in S-factories. 

A critical study of factorization assumption is an urgent task, and has indeed 
been started ||, at a time when increasingly accurate experimental results learn us 
that our present understanding of B — > ipK(K*), based on factorization assumption 
among other hypotheses, severely fails |4j]. 

It is usually claimed that the corrections to factorization are due to final state 
interaction (FSI). In a sense this statement is true, but the proper meaning of what 
is understood by that needs clarification. One of our aims in this paper is to proceed 
to this clarification in a simple model in which the dynamics is rather transparent. 

Furthermore little is known about the validity of factorization except that it is 
violated in low energy K — > tttt and D — > Kir, KK, txti |J |j| channels where strong 
FSI phases are experimentally known to be present, both from the direct analysis 
of weak nonleptonic decays and from scattering experiments. The study grows more 
difficult at larger energies when more channels are coupled. This happens when mul- 
tiparticle channels come in, particularly multi-pion channels, and also when several 
two-body channels communicate via strong interaction through quark exchange or 
quark-pair annihilation-creation. 

We will concentrate on two-body channels communicating via quark exchange. 
In a Tamm-Dancoff type expansion, quark exchange is dominant since it needs no 
quark pair creation neither annihilation. Furthermore, it has been stressed by J. 
Donoghue |7j that such a mechanism might explain the D° — > cf>K decay amplitude. 

On the other hand, M. Shifman || (see also [[J]) has made the interesting 
remark that beyond leading 1/N C order, the factorization assumption simply violates 
duality. As we shall show at length in this paper, this effect is fully related to final 
state interaction via quark exchange. 

Indeed this comes from the fact that in color space a qqqq color singlet can be 
decomposed in two ways into two qq color singlets, but the two resulting states are 
not orthogonal to each other ||. At the time the weak decay of the initial meson 
creates a qqqq color singlet, the four quarks interact strongly with one another. It is 
impossible to tell that one pair qq is in one meson and the other pair in another meson. 
However the system evolves, and eventually splits into two qq color singlets that are 
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spatially distant and hence orthogonal. During all the period of the interaction, one 
given quark does not know with which antiquark it is paired. This is exactly the 
situation which is expressed usually by the expression "quark exchange". What is 
depicted by this expression is not a simple and instantaneous exchange of a quark 
from one meson to another, but a period of overlap of the two mesons resulting finally 
in a non- vanishing amplitude for an exchange of a quark. 

The factorization assumption totally overlooks this complex mechanism, as 
it simply computes the overlap of the qqqq system resulting from weak decay with 
the final mesons. Two non- orthogonal states are taken as an approximation of two 
distinct final states which are obviously orthogonal. Once one realizes this, it is not a 
surprise that one encounters some problems with probability conservation. 

Our aim is mainly to understand better this interaction mechanism in a simple 
model, to check that when the dynamics is correctly treated there is no contradiction 
with duality, and to identify the different effects contributing at next to leading order 
in l/N c . 

We will work in the kinematical situation considered by Shifman ^, in a 
Hamiltonian approach, namely a quark model one. In section 2. we will rephrase 
duality in the language of closure theorem and reexpress the contradiction between 
factorization and duality. In section 3 we will present a quark model, with an adi- 
abatic approximation and mass-degenerate antiquarks, and we exploit the resulting 
permutation symmetry to compute the S matrix and the weak decay amplitude. In 
section 4 we will conclude. 

2 Duality versus factorization, rephrasing the prob- 
lem. 

Let us recall the ideal process which is studied in ||. There are three heavy quarks 
A, B, C and two light antiquarks a, f3. We will assume the latter to be heavy enough 
to justify the use of the quark potential model which will be our tool all along. 
The process under study is the weak decay 

P A5L -> Mbol + M c - p , M B0 + M CBl 

where Paq is pseudoscalar meson composed by A and a and M BBl represents any 
meson composed by B and a. 

The following relations are assumed ||: 

M B M r - M. M A = 2M + A, A QCD < A < M (1) 
To which we add 

Aqcd < m a , mp (2) 

to justify the use of quark model. 
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We cannot find a physical example of such a situation. If the s quark was 
heavy, the ideal situation assumed in |§ would be realized by the couple of decay 
channels: B d -> D°K° , D+R- . 



2.1 Relation between duality and closure theorem. 

In || M. Shifman exhibits a contradiction between duality and the standard factor- 
ization hypothesis also encountered in [[J while studying Ar for the B s — B s system. 
We will study this issue as a contradiction between the closure theorem and the stan- 
dard factorization hypothesis. Indeed, duality is related to the closure theorem in 
quantum mechanics as we will now recall. 

Let us call generically \n > all hadronic states built up with quarks B, C, a, j3. 
Calling Hw the weak Hamiltonian, the state Hw\Paq. > is composed of the four 
quarks B,C,a,/3. The decay width of the Pau meson is given by 



T(Paq) = En < PAa\H w \n >< n\H w \P ASl > 5(E n - E ASt ) 

= En < PAa\H w \n > 8(H - Eao) < n\H w \P A a > (3) 

where E A a is the initial energy, E n is the energy of the state \n >, H is the strong 
Hamiltonian (H\n >= E n \n >), and where the sum is to be understood as a sum 
over discrete states and an integral over continuum states. The set of states \n > 
is a complete set. The sum could be expanded on any basis, and in particular on 
the basis of the free quarks B,C,a,(3. This is where the closure of the Hilbert space 
comes in. 

Now, H = H c + V where H c stands for the kinetic energy and V for the 
potential. Whenever the contribution from V can be neglected in front of H c , one 
recovers lowest order duality, i.e. the simple parton model with no perturbative 
corrections neither non-perturbative ones from higher dimension operators. We do 
not want to go into the question of when this approximation is valid, and how to 
get better approximations. We simply want to rephrase the contradiction between 
factorization and duality as a contradiction between factorization and closure, then 
go into a simple model to show how the dynamics solves this. 

2.2 Contradiction between factorization and closure. 

Let us consider for example the weak Hamiltonian: 

H w = 2 V2G [d (B^LA) (C^L[3) + C 2 (B-y^P) (C^LA)] (4) 

where L — (1 — 7s)/2, C\ and C 2 are coefficients that we do not need to specify in 
this paper, although they are reminiscent of the familiar coefficients in the effective 
weak interaction Hamiltonian. 
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At the time t = the weak Hamiltonian acts on the initial meson and produces 
a state that contains the quarks C,/3,a,B. Let us call this state |/ >. We find it 
convenient for later use to decompose \f > into its color part and the remainder: 

|/ >= H w \P A a >= J dPcdppdp B dpo&(pc, sc,Pp, Sp,p B , S B ,Pa, So) 

$B, Sq, C B , Cq 

^[Cibcc^cBfia + C 2 5 CCtC J CBiC p} \C,p c , s c , c c ; (3,Pp, sp, cp, B,p B , s B , c B ; a,p & , s & , c a 

where pc, sc, cc (f>p, sp, cp) labels the momentum, spin and color of the (anti-)quark 
C (/?), N c is the number of colors. The function \l/ may be computed from the wave 
function of Pa q and the operator H\y in (Q). However we will skip this computation 
since the precise expression for \1/ is not relevant for our argument. 
It is obvious that: 

</!/>= (Cl + Cl + 2^) K (6) 

with 

K= Yl J dp c dppdp B dpa ^(pc,s c ,Pp,Sp,p B ,s B ,Pa,Sa) 2 (7) 

SC, Sp, C C , Cp, 
S B , $a, CB, Cg 

Let us also decompose the meson wave function into a color part and the 
remainder: 

dpcdppi^pipc, sc,Pp,Sp) — 5 Ca , C p \C, pc,sc,c c ;f3,Pp,Sp,Cp> 

(8) 

and analogously for all quark-antiquark pairs. The f/^j's form a complete orthonor- 
mal basis of the spin-momentum Cf3 Hilbert space. Let us now define the spin-space 
overlaps: 



K chB a = £ sc , s 0,s ,sa I dp c dppdp B dpa &(pc, S C ,Pp, Sp,p B , S B ,Pa, So) 

^cl(PC, Sc,Pp, S$)ll>M(fiB, S B ,Pa, S & ) (9) 

and analogously for the alternative grouping of quark-antiquark pairs: Ca; B(3. Clo- 
sure in the C/3 and Ba spin- momentum subspaces implies that 

Z^ £ ^Cf3;BaCI3;Ba 1 1 ' ^ ±u / 

n,m 
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and analogously 

n,m 

From eqs. (|]), @ and (Q) it results that 



< J\ M C^ M Ba>- a l K B frCa 



<f\M^;M^>= a 2 K^. (12) 



with 

ai = Cl + F' ° 2 = ° 2 + F (13) 

Up to now all equations were exact. Now we shall formulate in our formal- 
ism the factorization approximation by assuming that the decay amplitudes are well 
approximated by the overlaps: 



T (P A& -> M^-mS) ^< f\M$M%$ >* (14) 

More precisely, the standard factorization assumption Q uses eqs. ( |12|) (|13|) 
and (0) with Ci, C% computed from the electroweak theory complemented with QCD 
radiative corrections. Bauer, Stech and Wirbel [[| have proposed a phenomanological 
factorization assumption that keeps eqs. ( |T2| ) and (O) but with ai and a 2 fitted to 
all known D, respectively B, decays. 

Within standard factorization, summing over all two meson final states and 
using eqs. (§), ©, (0), © (|IJ and flTJ we get: 



E \T {Pa, - MgMg) | 2 +|T (p Aei - M^g) | 2 = ( (C? + C 2 )(l + -1) + 4^ ) X 

n,m \ ^ c c ) 

(15) 

To leading order in l/N c eq. ( |15f) gives the same result as eq. (||). In our 
present framework this reflects the well known fact that to leading order in 1/N C 
factorization and duality are compatible. However, the 0(1/N C ) corrections show a 
discrepancy, the contradiction stressed in g. As suggested in |§, this discrepancy 
could be cured, except for the 1/iV 2 terms by using a phenomenological factorization 
with 

- = 0. + ^, «i = C 2 + §- (16) 

However, in our framework it is easy to trace back the origin of the discrepancy. 
The fact is that the set of states \M^lM^ > @\M^IMq^ > is not an orthonormal 
basis of the Hilbert space. The states are normalised, but they are not orthogonal: 
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< M^M^W^M^ >= O(i-) ^ (17) 

in general. This overlap is 0(1/N C ) as can be easily derived from the color part in 
wave function (|8T): 



.Nc ) Cc,CB,Ca,Cn 



^Cc,c Sc B) c s Sc B) c 9 Sc c ,c a 
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It is easy to check that ( |TT| ) is at the origin of the discrepancy between (|T5| ) and (||). 

Neither are these states eigenstates of the strong Hamiltonian. Indeed, these 
states are built up from two asymptotic mesons combined via a plane wave for the 
relative momentum between the two mesons. When the two mesons lie far apart, 
the simple product of their wave function is an eigenstate of the strong Hamiltonian. 
However, in the states we consider there is a non negligible contribution with the two 
mesons overlapping in space where they strongly interact leading to an important 
distortion from the simple product of asymptotic meson wave functions. 



3 An Adiabatic Quark Model with Degenerate An- 
tiquarks. 

3.1 The S matrix 

Let us first consider an oversimplified model. We will assume all quarks to be spinless 
and a and (3 to be degenerate in mass: m s = mp = m. Next, B and C being very 
heavy, we will treat their motion as classical. They are supposed to move head-on 
with velocity v. As a function of time t the spatial coordinates of B and C are 

r B = -re = vt. (19) 

From now on the mesons will be assumed to have their center-of-mass localised 
in configuration space, at the position of the heavy quark and we will neglect 0(1/ M) 
corrections. The fact that the heavy quarks meet at the origin, i.e. that their impact 
parameter is zero means that our model will describe the S wave channel. 

Concerning the Hamiltonian for the antiquarks (3 and a we will use a color 



potential introduced in flO|| : 



H(t) = £ + f| + E ^bKV% - vt) + X a c X%V(f^ + vt) + \ a B \%V(r p - vt) 



A&ASV(fg + vt) + \%\%V(rs - r & ) + X a c X a B V(2vt) (20) 
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where \q is the Gell-Mann SU(3) Hermitean matrix applying to quark C, etc, and 
where — V(r) is a rotation-invariant confining potential, so that color singlet mesons 
are bound together (the A a A a factor is negative on a singlet). The Hamiltonian is 
bounded from below when restricted to overall color-singlet states. 

H(t) in is invariant for the Permutation : P = at <-> (3. It results that all 
eigenstates of H will be eigenstates of P with eigenvalue ±1. The asymptotic states, 
when T — > ±00, are built from simple products of the mesonic wave functions whose 
center of mass are located at ±vT: 



V2\D ± > n > m ,T> = \M^(-vT) > ®\M { j£{vT) > ±\M^(-vT) > ®\M ( ™\vT) > 

(21) 

where the M^ n \f) states are the mesons states defined from eq. (|8|) by a Fourier 
transform on the center of mass variable. We will assume the evolution in time to 
be adiabatic, i.e. we assume that the state \D ± ' n ' m ,t > evolves in time by remaining 
an eigenstate of H(t) for all t. We will further assume that during the evolution, 
the fundamental states \D ±,0 '°,t > never cross other states. It results that the two- 
dimensional subspace spanned by \D ±,0,0 ,t > is stable under the action of the strong 
Hamiltonian, i.e. that any state within this subspace evolves into a state within this 
subspace at a later time. Consequently, the two-by-two restriction of the S-matrix to 
this subsapce has to be unitary. We will now restrict ourselves to the study of the 
strong interaction scattering process of these two fundamental states. At time t, the 
eigenstates verify: 

H^D^^t >= E ± (vt)\D ± ^°,t > (22) 

where we made use of the fact that H depends on t only through the product vt. The 
interaction between the two terms in the r.h.s. of (|2lD is 0(1/N C ) as already argued, 
see flTB|). Hence 

E + (vt) -E~{vt) =0{^-) (23) 



•N, 



c 



Asymptotically it is also obvious that 

E + (±oo) = £T(±oo). (24) 



since the two mesons do not overlap implying that the two terms in the r.h.s. of (p2T|) 
become orthogonal. 

In fact the two energies differ from zero only when the two mesons overlap and 
the overlap falls off exponentially when vt — ► 00. In the basis \D ±,0 '°, ±00 >, (plD , 
the S matrix is diagonal and, being unitary its general form writes: 

S = e2tS ( e 7 e-«*) (25) 
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where 



25 



+°° dz 

oo V 



E + (z) + E-(z) E + (oo) + E~(oo) 



2 2 

dz E + (z)-E-(z) 



-oo V 

Indeed, the S-matrix is given by 



S = Texp 



dtH^t) 



(26) 



(27) 



which in our case of a state, say |+ >, that remains eigenstate of a time dependent 
Hamiltonian with energy E + (t), simplifies to 



exp 



/oo 
dt(E + (t) - E + (oo) 
-oo 



(28) 



where the interaction Hamiltonian has been taken to be the total Hamiltonian minus 
the Energy of two non-interacting mesons. 
In the meson-meson basis: 



the S matrix writes: 



m^I{-vt) > ®\Mzm > mi{-vT) > ®\M%m > 



D 2i8 ( cos 20 % sin 20 
I i sin 20 cos 20 



(29) 



(30) 



This matrix is, as expected, unitary and invariant for the permutation P which 
permutes the lines and the columns of the matrix. When the angle does not vanish, 
there is a quark exchange between the mesons which becomes maximal for = tt/A. 
Actually, = 0(1/ N c ) from eq. (p3|) . From eq. (|26D , we see that oc 1/v. The 
reason for this is clear: the lower the velocity, the longer the mesons overlap and can 
exchange quarks. We will now return in our model to the contradiction discussed in 
the preceeding section between duality and factorization. We work out an illustrative 
example of these features in the Appendix. 



3.2 Final state interaction in our model. 

In this section we return to the weak interaction. We consider here the exclusive 
decay channels Paq — * M^M^l and Pao — * Mq\m^\. Thanks to the statement 
made in the preceeding section that the subspace spanned by these two final states is 
stable for the weak interactions, we can safely forget all other channels in our study 
of the FSI. 

We now assume that the weak Hamiltonian acts at t = 0, creating the heavy 
quarks B and C at f — 0. The weak interaction creates a state |/ > as defined 
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in (|jD. When t = an additional symmetry is present in the strong Hamiltonian 
H(0): invariance under permutation of color labels Cc Cb and consequently 
under Cp <-> C s . It results that the |.D ±,0 '°,0 > states are even (odd) under the 
color permutation P c : Cc ^ Cb- Restricted to the color-even (color-odd) sector the 
Hamiltonian H(0) reduces to: 

where the upper (lower) sign corresponds to color-even (color-odd) states. 

It is important to notice that to leading order in N c the Hamiltonian (|3l|) is the 
same for color-even and color-odd states. It results that in the iV c — > oo limit, the two 
color wave funtion multiply the same spatial wave function for t — 0. Furthermore, 
to leading order in N c the Hamiltonian fl31|) is equal to twice the Hamiltonian for one 
heavy-light meson: 

p L _2 W -l) W+ l) 

2m iV c v ; v ' 

where p and f are the light quark momentum and position. This means that to 
leading order in N c the Hamiltonian ( |5T] ) corresponds just to the sum of two non 
interacting mesons superposed at the origin. This corresponds to the factorization 
assumption. 



Eq. (|3T| ) exhibits a symmetry for the exchange of spatial variables 



r 



0- 



This symmetry is simply a product of the color permutation symmetry P c , valid at 



t — 0, and the global permutation symmetry P, valid for all t. The eigenstates of ( 51 ) 
are eigenstates of the spatial permutation r a — > rs, and it is not difficult to guess that 
the ground states are symmetric under the latter permutation. This is illustrated in 
the Appendix. 

Hence, restricting ourselves to the subspace spanned by the two fundamental 
states, which are symmetric states for the permutation r*„ — > fa, we have P c = P and 
we project the state \f >, ©, into the subspace 

Ho = ID+' ' , > ®\D-'°'°, > (33) 

We have: 

< D^Olf >= v ~^~; /v 2 ~ x 7 — S ± (34) 



where S 11 * 1 is the spatial overlap. As stated above S ± = S° + 0(1/ N c ) where 5*° is the 
spatial overlap of \f > with the direct product of two non interacting mesons located 
at the origin, i.e. ip^ (fa)^ 1 ' (fa) , ip^(f) being the ground-state eigenfunction of 

The evolution forward in time of the states \D ±,0, °, t > is obtained by replacing 
in ( p6|) the — oo lower bounds of the integrals by 0. It results, thanks to time reversal, 
in phase shifts which are simply divided by 2: e tS±t ^. 
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The resulting T matrix for the decay of the initial meson Pao. into the funda- 
mental mesons is: 

T{P A& - M^MSl) = e» | (Cl+C2)(1+ t )S+ e^ + ^-^-fo 5 ' e -*\ 

Cl3 1 2(1+^)2 2(1-^)3 J 

To perform a systematic l/iV c expansion, let us first define AS 11 * 1 by 

S± = S° + ^ (36) 
Then, from (|35| ) we obtain to first order in 1/N C : 



Cf3 lv± Ba) ^ >-> \ \y± ~ 2N C ' 2N C S° 

+H^ + ft + g2(A 2S AS " ) ) Si ^> 

T(P^ - MgjM^) = C «S°{ (<7 2 + + ^±^) cos 

+^fci + ^ + Cl(A 2 5 ; e 5^ ) )^^ (37) 

Comparing ( p7|) with eqs (|T2|) and (|14D we see that if we take in (^) AS" 1 * 1 = 
= we recover the factorization corrected a la Shifman, i.e. with eq. (0). The 
global phase 5 is not relevant here. Consequently we learn that this factor 1/2 a la 
Shifman is completed by two other effects at the same order in 1/N C : i) the AS"* 1 
which reflect the difference, at the time of the weak decay, between the total qqqq 
spatial wave function and the simple product of the two asymptotic meson wave 
functions; ii) the phase which reflects strictly speaking the final state interaction. 
The latter phase 0, although 1/N C suppressed may become very large when v — > 0. 
For = 7r/4 the factorization becomes grossly wrong since the role of the two final 
states is interchanged: The operator multiplying C\ in Hw produces dominantly the 
M^Im^ instead of M^M^l as suggested by factorization, and vice versa for C^. 
In other words, the class II decays become dominant over the class I when = 7r/4. 

Furthermore, the dynamical origin of the phase and of the factors AS' ± is 
obviously strongly dependent of the precise nature of the decay channels considered. 
This indicates that the effects of and AS^ cannot be incorporated in a phenomeno- 
logical factorization a la BSW, which assumes a given pair of constants a± and a-i for 
all the decay channels of a meson. 

In the appendix we have performed an explicit calculation of AS^ and an 
estimation of for an harmonic oscillator potential. As can be seen from ([71]) and 
|73| ), is always the dominant 1/N C correction when a and (3 are light quarks, and 
when they are heavy, still dominates as long as vmR < 1. 

Let us now consider a more general case by first introducing spin. 
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3.3 Final state interactions of the fundamental pseudoscalar 
and vector mesons. 

From heavy quark symmetry (HQS) we know that vector and pseudoscalar mesons 
are degenerate. HQS also tells us that the heavy quark spin is conserved. This is 
of course a trivial consequence of fl2"0|), but it is quite general. For example, in (p0|) 
we might add a term coupling the light quark spins, such as ■ <Ja, but all terms 
including the spins sb and sc are 1/M suppressed. 

Restricting ourselves to a total quark spin, the ground states combine into 
four possible asymptotic states: P C gP Bdl , Vc^Ba-, Pb^Pco and V B nVca (where P 
stands for pseudoscalar and V for vector). 

It is then convenient to use states with a given symmetry for Pg = sb <-> sc- 
The relevant combinations are: 

Ps { 2 1 = "1 — » — 



where W stands for V°V'° — V + V'~ — V~V' + with 0, +, — labelling the polarization 
of the vector mesons. 

In fact, the first combination in (|38|) corresponds to Sbc — (total spin of B 
and C) and the second to Sbc — 1- Using for large |T| the notation: 

ii n (o,o) T f -|P C/ j(-gT)P aa (^)>+|Vb / j(-^r)VB a (^)> 1 

I C/3;Pa' 1 2 J 

\ 3 D {0 '°KT>= I - 3 l P C/3("^) P aa(^)>-|^Cg(-^)^(^)> 1 (gg) 

the strong Hamiltonian is diagonal in the basis where both Ps and P = a ^ (3 are 
diagonal: 

vW°,T> ^ | 1 B« a ,r>±|'D« J ,T> 



The four states in fl40|) evolve diagonally under the strong Hamiltonian and 
lead to four phase shifts given by formulae similar to eq. fl2T5p . 

Next we make the assumption that \f > defined in (||) is odd under Pg: 

P s \f>=-\f> (41) 

Relation ( flip is a consequence of Fierz symmetry whenever Hyy is build up 
of Fierz-invariant currents, as is the case in eq. (HI). The fact that Fierz symmetry 
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translates into a spin antisymmetry as in ( f4T|) comes from the fact that Fierz transfor- 
mation contains an additional minus sign from fermion field commutation. It results 
that only the states I 1 /} 1 *"' 0,0 , > are produced during weak decay, i.e. Sbc = 0. The 
arguments from the beginning of section [372] to eq. ([34]) may be repeated, except that 
due to the spin asymetry, the \ 1 D + ' 0,0 ,0 > (l 1 /) - ' 0,0 ^ > is color-odd (color-even): 



< D ' ' 0\f >= — — 2 a ^- ST (42) 



leading to 



1 CI3 ' BaJ \ 2(l-£)* 2 (l+£)* / 

V Cq;B/3; \ 2(i-^)4 2(i+£)i J v ; 

and to first order in l/iV c : 

7YP< _^M (0) M (0 h- VpiSoOffr _l_ C 2 , Ci(AS++A5-) 



L C/3 ±v± BaJ 2^ " I V^ 1 ' 2Af c 1 2iV o 5 y COS ( 

-,(C 2 + ^ + ^^-) )sin0} 
T(P AS - MSJM^) = fe^°{ - (C a + ft + cos0 

^fci + ft + Cl(A 2 5 ; c ^' ) )sin0} (44) 

where r\ = +1 for longitudinal vector mesons and rj = — 1 for transverse vector and 
pseudoscalar mesons. The difference between the r.h.s of (H31) and the r.h.s. of 



comes from the interchange of color symmetric with color antisymmetric combinations 
as apparent when comparing (|34]) and (f£|). 

It is to be noted that the relation between PP and VV production ampli- 
tude is exactly given by the fact that only the X D combination are created. Indeed, 
this relation is a consequence of HQS and would only be corrected if we considered 
the 0(1 /M) corrections. This relation is not a surprise when one realizes that the 
conditions ([!]) and (0) imply an S-wave dominated decay. 



3.4 Semi- inclusive decay. 

We have found important channel dependent corrections to factorization. We may 
still wonder if these corrections ares not washed out when we sum up on one side 
all the decay channels M^M^ for all m, n, and on the other side all the channels 

m bI m S!- This is the 

aim of this section. 

Let us call TC + (t) (H~(t)) the Hilbert space spanned by the set of states 
\D+' n ' m ,t >,Vn,m (\D-' n ' m ,t >,Vn,m) defined in eq. (P). H + (t) contains 
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the even (odd) states under the permutation P = a «-> (3. The latter commuting with 
the Hamiltonian H(t), the evolution does not mix the spaces Ti. + (t) and H~(t). We 
shall call U + (jt 1 ,t 2 ) (Z7~(*i, t 2 )) the evolution operator in H + (t) (H~(t)). U ± (t 1 ,t 2 ) 
are unitary. 

As already stated in section |3T^ , H(0) is also invariant under the color permuta- 
tion P c : Cc <-> Cb. For £ = the Permutation P = P c P r where P r = p & ,s & ps, Sp 
(remember r B = r c = for t = 0). We then decompose \f > both into eigenstates 
of P: 



I r+ ^_ (Ci+C2)(i5c C ,c^c B ,^+^,^^,^1 , s (Ci-C2)(<5c c ,c^<5 Cs , C Q-<5 cc , C(i (5c s , c ^) . s 

I »_ ( c i- c 2)(Sc C ,c^Sc B ,o a -&c C , aBl 8c B ,cp) i 5 (Ci+C2)(<5 cc , c ^<5 cs , C[i +(5 cc , C(i 5 cs , c ^) „ 

I-/ >— 2N C \J+ > H 2iVc \J- > 

(45) 

where |/ ± > are eigenstates of P with eigenvalue ±, \fs> contains the spin-space 
part of the wave function (|5|): 

l/ S >= £ _ _ Jdpcdp0dpBdp & V(pc,Sc,P0,Sp,PB,SB,Pa,S a ) 

s Ci Sp, Cc, Cp, r 

0b> Ca 

\C,Pc,SC,Cc]j3,Pp,Sp,Cp]B,pB,SB,C B ;a,Pa,Sa,Ca > (46) 

which is expanded into the eigenvectors of P r , \f± > corresponding to eigenvalues 

P r \fi>=±\fi>. 

The norm oi \f^ > is: 



</ ± l/ ± >4 



{Ci ± so 2 (i ± y) < fK > +(^ * oo a (itj)< > 

(47)" 

If we define < /f |/f >= if±, then K + + K_ = if as defined in eq. (0). 
As stated above: 

I/* >e (48) 
The evolution toward a large positive time T leads to: 

\f ± (T)>=U ± (T,0)\f ± > (49) 

By unitarity the norm of \f ± {T) > equals that of |/ >. 
We get 

V2T(P A& -> D±> m ' n ) = (d ± C 2 )(l ± S^ m ' n e^+ m ' n + 

(d =F C 2 )(l =F f )5JTE S^V*^" (50) 
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where S±' m ' n and 0jjj' m ' ra are real numbers denned by: 

1 



< £) ± > m > ri ) T|/^(T) > 



with 



£(T) = 1^(^0)1^(0) > 

(/_) being the first (second) terms in the right hand sides of (f45|). 
From unitarity: 



which leads to 



±,m,n 



E \T(PA, - D^)\ 2 = (Cl + Cl + 2^) 



as expected from @ and unitarity. 
Finally 



(51) 



(52) 



(53) 



(54) 



T(P 



As 



T{P J 



Aa 



K 2 
K 2 



\ (C 1 +C 2 )(l+^-)S 



1 \ c~l~' m ' n 



i \ ( Cl -c 2 ){i-^- c )st' m ' n 

i f (Ci+C 2 )(l +1 i-)S+' ra '" 



+ 
+ 



K 



\ \ (C^C 2 ){l-^- c )St' m ' n 



(Ci 


-c 2 )(i-^)5;^ 






(Ci- 


fc 2 )(i +1 ^)5:- m -" 






(Ci 








(Ci- 


fc 2 )(i +1 ^)5:- m '" 


2(i +1 ^)i 



(55) 



Duality is fully verified since summing over all states 
E \t{Pa* - M^Mg)| 2 + |T(P Aa - MgMg)f : 

let us now consider the partially inclusive sums: 



iV, 



(56) 



E \ t (Pa* 

m,n 

E \ t (Pa^ 



M^Mtl 



= £ 



C/3;Ba 



J Ca;Bf3 



(57) 
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Nothing general can be said. In the iV c — > oo limit, 

n+,m,n n—,m,n ±+,m,n i—,m,n /ro\ 

S± ^ S± ; 0±* ~ 0±' (58) 

Assuming that for finite N c we keep the relations (|58|) , we recover Shifman's 
Ansatz: 

= f {(Ci - c 2 )(i - + (d + c a )(i + ^)^} 2 = if {cl + + o(^)} 
Wis = f {(Ci - c 2 )(i - - (Ci + c 2 )(i + = if {c 2 2 + + O(^)} 

(59) 

If we made the opposite assumption that the relative phases 0±' m,n — 0±' m,n 
are random, which might be reasonable at small velocity, the result would be 

E C ft* a = = * { + ^ 1 (60) 



that is, an equal sharing of the total probability between the two sets of channels. In 
the latter case, even though we consider an inclusive sum, the final state interaction 
has a non trivial effect: the quarks have been redistributed at random between the 
final mesons. Of course such a random phase equal sharing may only happen when 
phase space allows for many final states to add up in a random way. 

As a side remark we would like to mention another, not yet published, study 
that we have performed on duality versus factorization. We have considered a model 
with non-relativistic scalar quarks bound to color singlets by a color harmonic oscilla- 
tor potential without assuming heavy quarks neither using an adiabatic approx- 
imation as done in the present paper. This model also automatically restores duality 
i.e. the conservation of probability. Summing over all mesons in the limit in which 
the radius R — > oo one should find the free quark result. This is indeed the case: 

EcftBa oc h + ^ A r] 2 + °( A ^) ( 61 ) 

2Wi9 oc [C 2 + (1 - y)^-J + O (J^\ (62) 

with y depending non-trivially on the masses, unlike the universal factor 1/2 proposed 
by Shifman. It is seen that FSI restores automatically duality. Notice by the way that 
expression (|62"D , for an arbitrary value of y, is more general than Shifman's Ansatz 
and would as well restore duality. 

4 Conclusions. 

We have used a Quark Model where the motion of heavy quarks is treated as classical 
and where we assume two mass-degenerate antiquarks. We have used the resulting 
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permutation symmetry to simplify the problem. We have restricted ourselves to the 
qqqq sector, and we have shown that the contradiction between standard factorization 
and duality stems from the non-orthogonality in color space of the two decomposition 
of the qqqq singlet into two pairs of qq color singlets. Taking care to use an orthogo- 
nal basis that diagonalizes the Hamiltonian, the dynamics of this sector shows very 
clearly how the final state interaction corrects standard factorization such as to satisfy 
duality. 

Shifman |1 has proposed to correct factorization by a replacement of 1/N C 
by 1/ (2N C ) while keeping the phenomenological factorization a la BSW. We have 
shown that this effect is indeed present. However we find two additional effects to 
same order in 1/N C . One is related to the spatial distorsion of the meson wave 
functions at the time of the weak decay: the two mesons overlap in space and hence 
interact strongly. This has been expressed by our parameter AS^. The second and 
more important additional effect is the phase difference <fi between the permutation- 
even and permutation-odd states when evolving after the decay until the mesons 
are spatially distant. The latter effect is 0(1/N C ) but also 0(l/v) where v is the 
final meson velocity in the total rest frame, and also 0(1/ mR) where m is the light 
antiquarks constituant mass and R is the wave function radius. This phase shift effect 
should dominate in the small velocity regime for light antiquarks. For large velocity 
and for heavy antiquarks, it vanishes. 

We have seen that in the exclusive case, restricting ourselves to the ground state 
mesons, the phase shift (f) is the dominant 1/N C correction for small velocity and light 
antiquarks, and it may produce a total modification of the factorization assumption, 
which could, for <fi = it/ 4, be large enough to totally interchange the amplitudes of 
the two channels, and lead to a dominant class II decay. In an illustrative example 
treated in the appendix the t = wave function distortion , AS" 1 * 1 turns out to be 
small. We ignore if this is a general feature. Would it be so, it would indicate a 
validity of Shifman's Ansatz for large velocity and rather heavy antiquarks. 

In the semi-inclusive case we compare the total decay probability into two sets 
of channels that correspond to the two possible pairing of qqqq into two qq. Again, 
the small phase shift, small distortion limit amounts to Shifman's Ansatz, while the 
opposite, random phase shift limit, amounts to an equal sharing of inclusive decay 
probability into the two sets of final channels. Again the latter situation may be 
reasonable in the small velocity case provided many final states are kinematically 
allowed. 
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A An illustrative example. 



As an illustration of the section |3.2| , let us take for the potential V in fl20|) an harmonic 
oscilator potential: 

^ = - 4 W -W * (63) 
so that the ground state solution of ( p2[) is: 

^W(r) = —^e'^ (64) 

where r = \r\. 

For the spatial part of the wave function \f > we take in configuration space: 

Wo^?B,?a) = G^(fa)5 3 (f c )d 3 (f B )5 3 (f p ) (65) 

which expresses the fact that the weak operator is localQ and that the quark a is a 
spectator coming from the Pa& meson, i.e. in the ground state wave function. G is 
proportional to the Fermi constant. 

To compute the spatial overlaps S^ 1 in we need to know the ground state 



solutions of the Hamiltonian ([H]) with fl63|) for V. Let us change variables in (|31): 

1 



fg = R- - f, fa = R+ \r, 



Pp= -jP- P, Pa = \P + P, (66) 



leading to 



H(Q) = ^ + ^ + (N c ±2)(N cT l) f32 (iV e ±4)(iV eT l) 
y > m Am A{N*-l)mR± (iV c 2 - l)mR* ' K ' 

The ground state solution, i.e. the spatial wave functions of |Z} ± ' 0,0 ,0 > is: 

3 3 (iV e ±2)^(fa+rg) 2 (jV e ±4)5(r^-rg) 2 

^'"'"(ra, rg) = ( c )5( f c ^ ^ e *C*c±i)*rf e (68) 

(iV c ± l)4i? 3 7T2 

The overlap is given by 

_ £ 4§(AT e ±2)l(Ar e ±4)l ^gg^ 

(A r c ±2)3+(7V c ±4)2-+2(AT c ±l)iZ 2 ,R§7rl 



2 Remember that we assume here spinless quarks. 
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where we have in (|65| ) left aside the S functions related to the heavy quarks, since the 
latter are treated classically. One has also for the factorization hypothesis result: 



S° = G— 3— 3- (70) 

-K27T4 

which is the ^(^(O) coming from the overlap of S(fp) with , 1 (f^). The overlap of 
ip^ifa) with itself gives obviously 1. S° in ( [70]) is obviously the N c — > 00 limite of 
S ± in f|). 

To next to leading order in l/N c the calculation of AS' ± defined in (^) is now 
straightforward from (|69|) and (|70|): 

which turns out to be rather small, mainly because the normalization factor compen- 
sates for a large part the modification of the integral. 
From (|BTD we also get the ground state energy: 

^0)~<l^Y + (^) } (72) 



leading to 



2mR 2 \N c ±lJ \N C ±1 




B+( °)-^( 0) = ^ + °l^l < 73) 



To compute <j) we need, (|2q), to know E + {z) — E~{z) for all z 7^ 0. This is not 
so easy to compute. We will simply use ([F3|) to make an order of magnitude estimate. 
We will assume E + (z) — E~(z) to be equal by -E' + (0) — -£^(0) as long as the hadrons 
overlap, i.e. for \z\ < cR where c is some number of order 10. Then 

' 3C (74) 



mvRN, 



c 



Although purely indicative, this result, besides confirming that <fi oc l/vN c also 
learns us that mR is the dimensionless number that gives the scale. When a and 
/3 are light quark, it is known that mR ~ 1. Comparing ([74]) with (|7TD and to the 
1/ (2N C ) correction a la Shifman, we check that in the exclusive case, the phase shift 
is the dominant 1/N C contribution. 



The wave function radius for the wave function Q is < r 2 >= 3/2R 2 . Hence one might think 
of taking c ~ 3/2. 
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